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O ■ Abstract: 



Effective theories are helpful tools to gain an intuitive insight into phenomena governed 
by complex laws. In this work we show by means of Monte Carlo simulations that Z(2) 
spin models with only spin-spin interactions approximate rather well the critical behavior 
of £77(2) gauge theory at finite temperature. These models are then used to formulate 
■ an effective description of Polyakov loop percolation, which is shown to reproduce the 

CN . thermal critical behaviour of SU (2) gauge theory. 

^> 

oo ■ 1 Introduction 

o , 

The global Z(N) center symmetry plays a key role in the confinement-deconfinement 
transition of SU(N) gauge theories. The spontaneous breaking of this symmetry is re- 
sponsible for the phase transition, and the corresponding order parameter is the Polyakov 
loop [||, 0. In full QCD the dynamical quarks act as a (small) external magnetic field 
and the conventional formalism of spontaneous symmetry breaking does not apply. There- 
fore it seems helpful to consider an alternative approach to critical behavior of SU(N) 
^ | gauge theory, which may be more readily generalizable to full QCD. 

In ||] it was suggested that the deconfinement phase transition in SU(N) gauge theory 
can be viewed as a percolation transition of suitably defined Polyakov loop clusters. In 
0, such an approach has been applied to (2+1)- and (3+l)-dimensional ££7(2) pure 
gauge theory; the clusters are formed by joining nearest-neighbour like-sign Polyakov loops 
with a bond probability similar to the one used to define the "physical droplets" in the 
Ising model @ . In |5], || , the bond probability was determined from an effective theory of 
Polyakov loops, which was derived by a strong coupling expansion [[| . Hence the validity 
of this approach is limited to the strong coupling limit of SU(2), and the calculations 
were thus done on lattices with a temporal extent of only two lattice spacings. It is thus 
not clear how to formulate a percolation picture for £77(2) in the continuum limit. 

The aim of the present paper is to clarify whether the deconfinement phase transition 
in (3+l)-dimensional SU(2) gauge theory can be described more generally as a percola- 
tion transition. We will derive an effective theory for Z(2) spin variables containing up 
to 19 spin-spin couplings which are obtained by solving Schwinger-Dyson equations, as 
originally proposed in @, [J. This method is independent of the number N T of lattice 
spacings in the time direction and we will apply it to two different cases, N T = 2 and 4. 
The particular effective theory studied here is then shown to allow an equivalent perco- 
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lation formulation of the deconfmement transition. The corresponding clusters are, with 
good approximation, the "physical droplets" of SU(2). 

2 Percolation and the Effective Theory 

We are interested in constructing an effective theory of Z(2) spin variables {s n = ±1}, 
defined as the signs of Polyakov loops at the spatial sites n. The effective Hamiltonian 
Tl(s,ft) of the signs {s n } of the Polyakov loop configurations can be defined through the 
equation || || 

exp[H(s,ft)] = f [dU]Y[ 8[s n ,sgn(L n )} exp(-S w (ft)), (1) 

n 

where L n is the value of the Polyakov loop at the spatial point n and Sw the Wilson 
lattice action which depends on the gauge coupling ft. Eq. (P shows that all degrees 
of freedom of the original SU(2) field configurations are integrated out, leaving only the 
distribution of the corresponding Ising-projected configuration. 

The problem is now how to determine the explicit form Ti,(s,ft), starting from the 
original Polyakov loop configurations. In general, TC(s,ft) will contain infinitely many 
operators. However, in the vicinity of the critical point where the correlations become 
long ranged one expects that H(s, ft) can be well approximated by a finite sum of different 
terms [§]. So we write 

n(s,ft) = EE^Wn, (2) 

n a 

where K a (ft) are the couplings, 0„ are the spin-spin operators at site n and a specifies 
the specific spin-spin interaction; e.g., in the simplest case there would be only a single 
term KiO^ = Ki Y^t=i s n s n+(i, with \x denoting a unit vector in /x-direction. In contrast to 
0, we here do not include multi-spin operators in Eq. (§), because it is not known how 
to define clusters in this case. 

To calculate the couplings, we use the following set of Schwinger-Dyson equations 

(01) = -(02ex P (-2E^)), (3) 

a 

which are derived by exploiting the Z(2) symmetry of TC(s,ft) ||. Here < ... > denotes 
averages taken in the (3 + l)-dimensional SU(2) gauge theory. Eqs. establish a relation 
between thermal averages of the operators and the couplings « 7 . These equations are, 
however, implicit in the couplings. They can be solved by means of the Newton method, 
based on successive approximations ||. One starts by making a guess about the values of 
the couplings and derives from Eqs. @ another set of values for /c 7 's. After a sufficient 
number of iterations, the series of partial values for the k 7 's will converge to the solution 
of Eqs. d). 

We notice that the general set of equations (|3]) refers to a single point n of the spatial 
volume. However, the thermal averages are independent of the particular point n, so 
it doesn't matter where we decide to take the averages. We thus explore translational 
invariance and the Eqs. (3) should thus be understood as being averaged over the entire 
lattice. This reduces considerably the effect of thermal fluctuations and, consequently, 
the errors on the final k~'s. 
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With this we have all necessary tools to derive an effective theory for 577(2) from the 
Polyakov loop configurations. 

In general, the approximation improves the more operators we include in (|2|). The 
fact that one must restrict the choice to some subset of operators involves a truncation 
error in addition to the statistical one. The error due to truncation is difficult to control. 
It can, however, be estimated by calculating the parameters of the effective theory for 
different sets of operators. We also need to establish a criterium to judge how well 
the effective theory approximates the original one. A good option could be to compare 
average values we get from the configurations produced by simulating the effective theory 
with the corresponding quantities measured on the original Ising-projected Polyakov loop 
configurations. In particular, we used the lattice average of the magnetization m, 

m = ^fI5>;I, ( 4 ) 

i 

(V is the spatial lattice volume) for such a test. 

We point out that the approach we have described can be used for any value of the 
number N T of lattice spacings in the temperature direction, although we have applied the 
method here only to SU(2) in (3 + 1) dimensions for N T = 2 and 4. 

The magnetization transition in the Ising model can be described as percolation tran- 
sition of appropriately defined site-bond clusters @. The clusters are defined by the 
condition that two nearest neighboring spins of the same sign belong to a cluster with 
bond probability (or bond weight) 

p = 1 — exp(— 2k). (5) 

Here k is the Ising coupling divided by the temperature. 

It was recently proved that this result can be extended to other spin models [[□], [TD|, |T2 



In particular, if a theory contains only ferromagnetic spin-spin interactions, the above 
percolation picture can be trivially extended by introducing a bond between each pair of 
interacting spins, and a relative bond probability 

p a = 1 - exp(-2« a ), (6) 

where K a is the coupling associated to the a-th interaction of the theory. 

As we will see in the next section, all couplings of the effective theory indeed turn 
out to be positive. To carry out the percolation studies we will build the clusters in the 
Polyakov loop configurations using bond probabilities given by Eq. @, with K a values 
taken from the effective theory. Once we have grouped all spins into clusters, we calculate 



the values of the following variables [13]: 



The percolation strength P, defined as the probability that a randomly chosen lattice 
site belongs to the percolating cluster. P is the order parameter of the percolation 
transition. 

The average cluster size S, defined as 

Here, n s is the number of clusters of size s per lattice site divided by the lattice 
volume, and the sums exclude the percolating cluster. 
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To determine the critical point of the percolation transition, we used the method 
suggested in [[TJJ] . For a given lattice size and a value of (3 we counted how many times we 
found a percolating cluster. This number, divided by the total number of configurations 
at that (3, is the percolation cumulant. This variable is a scaling function, analogous to 
the Binder cumulant in continuous thermal phase transitions. Curves corresponding to 
different lattice sizes will thus cross at the critical point, apart from corrections to scaling. 

The percolation cumulant gives not only the critical point (3 C , but also the critical 
exponent v. In fact, if we rewrite the percolation cumulant as a function of (/3 — P^L 1 ^, 
where L is the linear spatial dimension of the lattice, we should get a universal scaling 
curve for all lattice sizes. We will adopt this method to determine the percolation exponent 
v. 



3 Numerical Results 

As we are interested in the phase transition of 577(2), we first focus on the derivation 
of the effective theory at the critical point. Simulations were done on 32 3 x N T lattices, 
with N T = 2 and 4. For the critical coupling we took the value (3 C = 1.8735(4) for 



N T = 2 as determined in §, and (3 C = 2.29895(10) for N T = 4 from fTJ. The Monte 
Carlo update we used consists in alternating one heat-bath and two overrelaxation steps. 
We measured our observables every 60 updates, which makes the analyzed configurations 
basically uncorrelated; the total number of measurements was 2000. Our aim is to check 
whether, at f3 — /3 C , we can find a projection of the theory onto the spin model defined 
byEq.©. 

For N T = 2 we considered two sets of operator containing 10 and 15 spin-spin couplings, 
respectively. The operators connect a point (000) to (100), (110), (111), (200), (210), 
(211), (220), (221), (222), (300), (310), (311), (320), (321), (322). In the case of N T = 4 
we have included further operators in our analysis, which connect the point (000) to (330), 
(331), (332) and (333). The final set of couplings is displayed in Table [I], where I refers to 
the projection with 10 (15) operators for N T = 2 (N T = 4) and II refers to the projection 
with 15 (19) operators for N T = 2 (N T = 4). Note that the physical length scale changes 
with N T . Two spins separated by a distance R measured in lattice units on the lattice 
with temporal extent N T are actually separated by distance r, with rT c = R/N T . Here 
T c = T(f3 c ) denotes the transition temperature of SU (2) gauge theory. 



From Table 1 one can see that the difference between couplings calculated for set I 
and set II is of the order of the statistical errors. This implies that possible truncation 
errors are small. Comparing our couplings for set I and iV r = 2 with the corresponding 
couplings obtained in || , we see that they agree within statistical errors, although multi- 
spin couplings were also included in the analysis of Ref. 0. This suggests that the 
multi-spin couplings are not important, for N T = 2; in fact, most of them are compatible 
with zero within statistical errors. Note that because of much improved statistics, the 
errors on the couplings in Table 1 are by an order of magnitude smaller than in |9j . Since 
the error on k 15 is of the order of its average value, we can set k 15 = without appreciable 
effects. In this way, we indeed have obtained the effective theory we were looking for, with 
only ferromagnetic spin-spin interactions. 
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Couplings 


iV T 


= 2 




= 4 




I 


II 


I 


II 


Kl (100) 


0.13042(8) 


0.13066(9) 


0.08385(3) 


0.08390(4) 


k 2 (110) 


0.01911(1) 


0.01905(3) 


0.01842(4) 


0.01839(5) 


^3 (HI) 


0.00476(5) 


0.00470(5) 


0.00769(5) 


0.00775(4) 


k 4 (200) 


0.00794(5) 


0.00801(11) 


0.00702(4) 


0.00697(1) 


^5 (210) 


0.00198(6) 


0.00192(4) 


0.00348(4) 


0.00343(2) 


^ (211) 


0.00069(3) 


0.00062(8) 


0.00203(5) 


0.00197(1) 


k 7 (220) 


0.00045(7) 


0.00033(2) 


0.00118(2) 


0.00114(1) 


k s (221) 


0.00013(6) 


0.00007(2) 


0.00082(1) 


0.00083(1) 


Kg (222) 


0.00017(4) 


0.00014(10) 


0.00038(7) 


0.00035(6) 


Kio (300) 


0.00072(16) 


0.00058(3) 


0.00113(4) 


0.00105(9) 


«u (310) 




0.00018(3) 


0.00076(5) 


0.00082(5) 


K12 (311) 




0.00008(1) 


0.00055(9) 


0.00055(4) 


km (320) 




0.00001(1) 


0.00032(3) 


0.00035(2) 


K M (321) 




0.00006(1) 


0.00033(1) 


0.00030(4) 


Kis (322) 




-0.00005(6) 


0.00029(3) 


0.00013(4) 


Kie (330) 








0.00020(5) 


«17 (331) 








0.00018(3) 


«is (332) 








0.00017(1) 


«ig (333) 








0.00017(4) 



Table 1: Couplings of the effective theory for the Ising-projected Polyakov loop configu- 
rations of (3 + l)-dimensional SU(2). I and II denote the two sets of operator containing 
10 and 15 operators for N T = 2, and 15 and 19 operators for N T = 4. 



Fig. [I] shows a comparison between the magnetization distribution for N T = 2 of the 
Polyakov loop configurations and the one of the effective theory containing 15 operators: 
the two histograms are very similar. The values of the average magnetization m are also in 
agreement: for SU(2), m = 0.091(1) and for the spin model, m = 0.0923(7). The average 
magnetization m calculated from the effective theory containing 10 operators is by 10% 
smaller than this value. One can also clearly see differences between the magnetization 
distributions calculated from the Polyakov loop configurations and from the effective 
theory with 10 operators. 

In the case of N T = 4 we have considered two sets of spin-spin operators containing 
15 and 19 couplings labeled I and II, respectively. The values of the couplings are shown 
in Table [l|. Note that all operators are positive in this case. The average value of the 
magnetization m obtained from the effective theory with 19 operators is 0.121(3), which 
should be compared with the corresponding value calculated from SU(2) configurations, 
m = 0.128(6). The effective theory with 15 operators gives a value which is by 20% 
smaller. 

The couplings reported in Table 1 decrease exponentially as function of rT. The 
exponential fall-off is approximately described by ~ exp(— 9rT) for both iV T = 2 and 
iV T = 4. It is easy to see that couplings of the effective spin theory at fixed rT should 
scale as N~ A close to the continuum limit. A comparison of couplings for N T = 2 and 
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N T = 4 at the same value of rT c shows that this scaling holds up to 20% scaling violation. 
This observation suggests that the form of the effective theory for obtained N T = 2 and 
4 will not change qualitatively by approaching the continuum limit. In particular, we 
expect that the effective theory will contain only ferromagnetic spin-spin coupling in the 
continuum limit. 

To check the volume dependence of the couplings we also determined them on 16 3 x 2 
and 16 3 x 4 lattices. The resulting values agree with those obtained from simulations on 
larger 32 3 x 2 and 32 3 x 4 lattices within statistical errors. We therefore conclude that the 
set of couplings we calculated is, with good approximation, also that of the effective theory 
in the infinite volume limit, which is what we will need for our subsequent percolation 
studies []. 

We have shown that the effective Z(2) theory with 15 (N T = 2) and 19 (N T = 4) 
spin-spin operators gives a good description of the thermal critical behaviour of SU{2) 
gauge theory. We now want to see if this effective theory also allows us to extend the 
Polyakov loop percolation description of the deconfinement transition 0, [5| from the 
strong coupling limit f| to a more general situation closer to the continuum limit. To 
study this we generate SU{2) Polyakov loop configuration on iV 3 x N T lattices, reduce 
them to Z(2) spin configurations and then define clusters with the bond weights as defined 
in Eq. (^. A cluster thus consists of a set of aligned spins coupled by one or more bond 
weights. We want to determine the percolation behaviour of these clusters. 

600 
500 
g 400 

c 
03 

t 300 
o 

O 200 
100 


-1 -0.5 0.5 1 

Magnetization 

Figure 1: Comparison of the magnetization histograms derived from the Polyakov loop 
configurations and from the effective theory (fj) defined by the couplings of Table |l|. 

We stress that the bond weights are temperature-dependent. Our effective theory 
represents a projection of 577(2) for (3 = (3 C . But in order to carry out our analysis, we 
need to evaluate the percolation variables at different values of (3. Strictly speaking, for 
each ^ we should derive the corresponding effective theory and use the set {^(A)} to 
calculate the bond weights (^) at /3j. We have calculated the couplings of the effective 
theory for several values of f3 in the interval [1.873,1.883] for 16 3 x 2 lattice. It turned 

x The bond weights are determined by the temperature. Once we fix the temperature, for different 
lattice sizes we must use the same set of values for the bond weights. 
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Figure 2: Percolation cumulants near the critical point for different spatial volume for 
N T = 2 (left) and for N T = 4 (right). The solid vertical lines indicate the value of the 
thermal threshold. For N T = 2 we also indicate the uncertainty on the threshold by two 
dashed lines. For N T = 4 the errors on the critical coupling are too small to be visible in 
the Figure. 

out that the variation of the couplings due to the variation of (3 in this small interval is 
of the order of the statistical errors. Because of that, at each f3, we shall use the same set 
of bond probabilities, determined by the couplings of Table [I]. 

For the analysis of the percolation transition we have used four different lattices of 
size N% x JV T , with N a = 24,30,40,50 and N T = 2,4. Fig. | shows the behavior of the 
percolation cumulant as a function of (3. In both cases the curves cross remarkably well 
at the same point, in excellent agreement with the thermal thresholds, indicated by the 
solid line. 

Next, we performed further simulations near criticality in order to study the critical 
exponents and to determine more precisely the value of the geometrical threshold. From 
standard finite size scaling fits of the average cluster size S at different (3 values we 
determine the value (3 C at which the scaling fit gives the best x 2 ; the slope of the straight 
line of the log- log plot of S at (3 C gives 7/V. The errors on f3 c and on the ratio of critical 
exponents are calculated by determining the /3-range containing (3 C such that for each 
(3 one still gets a good x 2 f° r the scaling fit. Unfortunately we could not determine 
P/u, because of strong fluctuations of the percolation strength P around (3 C . The value 
of P at criticality is in general quite small, and it is more strongly influenced by the 
approximations involved in our approach. Furthermore, we determine the exponent v 
from the finite size scaling of the pseudo-critical points, i. e. the positions of the peak of 
S for the different lattices we considered. 

The percolation threshold and the critical exponents are summarized in Table ^| for 
both N T = 2 and 4. There we also give the critical exponents of the thermal SU(2) 
transition and those corresponding to the 3D Ising model. As one can see from Table 2, 
good agreement between the percolation exponents and the 3D Ising exponents is found. 
In Fig. [3] we show the scaling of the percolation cumulant with the Ising exponent and 
with the 3-dimensional random percolation exponent v = 0.88 for the case of N T = 2. 
The figure shows that scaling is quite consistent with the Ising exponent while the random 
percolation exponent is ruled out by the data. The situation is the same for N T = 4 as 
one can see in Fig. fl. 
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Figure 3: Rescaling of the percolation cumulant curves shown in Fig. |2] for N T = 2 with 
the 3-dimensional Ising exponent v = 0.6294 (left) and with the 3-dimensional random 
percolation exponent v = 0.88 (right), using (3 C = 1.8734 for the critical coupling. 





Critical point 




V 


N T = 2 


N T = 4 


N T = 2 


N T = 4 


N T = 2 


N T = 4 


Percolation 


1.8734(2) 


2.2991(2) 


1.977(17) 


1.979(16) 


0.628(11) 


0.629(0.011) 


Thermal 


1.8735(4) 


2.29895(10) 


1.953(9) 


1.944(13) 


0.630(10) 


0.630(11) 


3D Ising 




1.967(7) 


0.6294(10) 



Table 2: Percolation critical indices for (3 + l)-dimensional SU(2). For comparison we 



also list the thermal results || and the 3D Ising values [16 . 



4 Conclusions 

We have shown that in the vicinity of the critical point SU(2) gauge theory can be well 
approximated by a simple Z(2) model with ferromagnetic spin-spin interactions alone. We 
have tested the accuracy of the projection for N T = 2 and 4. In both cases we have found 
that there is basically no trace of frustration. The percolation picture of the effective 
model is suitable to describe the deconfming transition of SU(2) in purely geometrical 
terms. 

In principle the present approach can be extended also to larger values of N T (smaller 
lattice spacing). We expect that for larger N T more operators should be included in the 
effective action. This can be avoided by first performing a block spin transformation of 
the Polyakov loop configurations and then doing the whole analysis for the blocked config- 
urations. We have shown that the couplings of the effective theory decrease exponentially 
in strength as function of rT c and their relative magnitude scales with N~ A . This sug- 
gests that the effective theory remains well defined in the continuum limit and will stay 
ferromagnetic and short ranged. Therefore we do not expect that the present picture will 
change qualitatively for larger iV r . 
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Figure 4: Rescaling of the percolation cumulant curves shown in Fig. |2] for N T = 4 with 
the 3-dimensional Ising exponent v = 0.6294 (left) and with the 3-dimensional random 
percolation exponent v = 0.88 (right), using (3 C = 2.2991 for the critical coupling. 
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